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Abstract 

We investigate the dynamic relaxation for SU(2) gauge theory at finite tempera- 
tures in (3 + 1) dimensions. Using the Hybrid Monte Carlo algorithm, we examine 
the time dependence of the system in the short-time regime. Starting from the or- 
dered state, the critical exponents 0, v and z are calculated from the power law 
behaviour of the Polyakov loop and the cumulant at or near the critical point. The 
results for the static exponents are in agreement with those obtained from simula- 
tions in equilibrium and those of the three-dimensional Ising model. The value for 
the dynamic critical exponent was determined with z = 2.0(1). 
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1 Introduction 



Traditionally, it was believed that universal scaling behaviour exists only in 
the long-time regime. However, recently Janssen, Schaub and Schmittmann 
[1] showed that far from equilibrium, in a macroscopic short-time regime of 
the dynamic evolution, there already emerges universal scaling behaviour in 
the O(iV) vector model. They considered the relaxation process of a system 
quenched from a disordered state to the critical point and evolving with dy- 
namics of model A (non-conserved order parameter and non-conserved energy 
[2]) and found that the magnetization undergoes an initial increase of the form 
M(t) ~ t e , where 9 is a new dynamic exponent and t denotes the time. This 
prediction was supported by a number of Monte Carlo (MC) investigations 
not only for O(iV) vector models, but also for several other systems with a 
second-order or a Kosterlitz-Thouless phase transition [3]. These simulations 
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offer also a possibility to determine the conventional (static and dynamic) 
critical exponents [4-6] as well as the critical point [6]. This may eliminate 
critical slowing down, since the measurements are performed in the early part 
of the evolution. 

In short-time critical dynamics exist two different time scales, the micro- 
scopic time scale i m i c and the macroscopic time scale t mac - Universal behaviour 
emerges only after t mic , the time during which the non-universal short wave 
behaviour is swept away. t mic is independent of the linear space extension L 
and it is in general small compared to the macroscopic time scale, which is 
proportional L z . After a time period of t mac the correlation length is of the 
order of L and the system crosses over to the long-time universal behaviour. 

First numerical simulations of the short-time dynamical relaxation at critical- 
ity started from a disordered initial state. However, dynamical scaling exists 
also for an ordered initial state. This is supported by a variety of MC investi- 
gations [3,7,8], but no analytical calculations exist for this situation. All static 
exponents as well as the critical point and the correlation length in the high 
temperature phase [9] can be also obtained with an ordered initial state. 

Systematic numerical simulations of the short-time critical dynamics have 
been carried out mainly in two-dimensional systems [3,10,11,8]. Also, a first 
approach to lattice gauge theory was done for SU(2) gauge fields in (2 + 1) 
dimensions [12]. The results strongly support that there exists a universal 
short-time scaling behaviour for the dynamic SU(2) lattice gauge theory. In 
this paper we investigate SU(2) gauge theory in (3 + 1) dimensions. In contrast 
to the previous study in (2+1) dimensions we start the dynamic evolution from 
an initial ordered state, i.e. with a magnetization of one at t — (mo = 1). 
For the determination of the critical exponents, the dynamic relaxation pro- 
cess starting from an ordered state has been proven advantageous over that 
from an unordered initial state. This is a consequence of less prominent fluc- 
tuations for an ordered initial state. However, the new critical exponent 9 can 
not been calculated if we start with m = 1. 

The dynamics of our system is given by the Hybrid Monte Carlo (HMC) 
algorithm [13]. Up to now mainly local algorithms such as the Metropolis or the 
heatbath algorithms have been used for the dynamic evolution. This is the first 
time that the HMC algorithm is studied in critical dynamic relaxation. The 
motivation for using the HMC algorithm is not only to enlarge the knowledge 
of short-time dynamics, but also the possibility to include fermions in future 
simulations. We examine the short-time behaviour of the order parameter 
and the cumulant at and in vicinity of the critical point. From the power 
law behaviour we extract the critical exponents (3, rj and z. The results are 
compared with those of simulations in equilibrium [14] and with simulations 
of the three-dimensional Ising model [7] , since it is expected that both models 
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are in the same universality class [15]. 

In the next section we sketch the scaling analysis of the short-time critical 
dynamics. Sections 3 and 4 briefly describe the model and the updating algo- 
rithm. Numerical results are presented in Sec. 5. The last section contains the 
summary and conclusions. 



2 Scaling relations 

Using renormalization group methods, Janssen et al. have shown that even 
in macroscopically early stages of a relaxation process 0(N) vector models 
display universal behaviour [1]. They studied a system initially in a disordered 
state with vanishing or small magnetization (m >0), suddenly quenched to 
the critical point and evolving with dynamics of model A, and derived the 
dynamic scaling form 

M W (t, t, L, mo) = b-W"MW (b~ z t, b^r, 6 _1 L, b x °m ) . (1) 

M(k) denotes the fcth moment of the magnetization, t is the MC time of the 
dynamic relaxation, r is the reduced coupling constant, b indicates a spatial 
rescaling factor and xq is a new independent exponent. For a sufficiently large 
lattice and small initial magnetization m t x °/ z this leads to 

M(t)~m t e , 9 = (x -f3/v)/z (2) 
at the critical point r = 0. 

Another important process is the dynamic relaxation from a completely or- 
dered state. For an initial magnetization exactly at its fixed point m — 1, a 
scaling form 

M {k) (t, r, L) = b' kfS/u M^ (b- z t, b 1/u r, b-'L) (3) 

is expected. The scaling form (3) looks the same as the dynamic scaling form 
in the long-time regime, however, it is now assumed already valid in the macro- 
scopic short-time regime. The validity of the scaling form (3) in the short-time 
regime was verified with MC simulations for a number different systems [3]. 

Taking b = t x l z for the spatial rescaling factor in Eq. (3) with k — 1 leads for 
the magnetization to a power law behaviour 

M(0~r ci , d = A (4) 
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at the critical point r = 0, if L is sufficiently large. For non-zero values of r, 
the power law behaviour will be modified by the scaling function M(l,t 1 ^ uz r). 
This can be used for a determination of the critical point. Also, the critical 
exponent \j{yz) [6] can be measured by taking the derivative with respect to 

T 

d T In M(t,r)\ T=0 ~t C;i , c,i = — , (5) 



while the dynamic critical exponent z can be determined from the behaviour 
of the cumulant 



Finite-size scaling shows that 



U(t)~t c v, c u = -, (7) 

z 



where d denotes the spatial dimension. Thus, the short-time behaviour of the 
dynamic relaxation starting from a completely ordered state is sufficient to 
determine all the critical exponents f3, v and z as well as the critical point. 
These measurements are usually better in quality than those starting from a 
disordered state. 



3 SU(2) gauge theory at finite temperatures 



The Wilson action for SU(2) gauge theory is given by 



where Up represents the usual plaquette term on the lattice. The number of 
lattice points in the space direction is L s and in the time direction L t . Thus, 
the volume and temperature are given by V = L s 3 L t and T — 1/L t , if we fix 
the lattice spacing a to unity. A point on the lattice has integer coordinates 
x = (x , x) = (xq, xi, x 2 , x 3 ), which are in the range < x < L t , < Xi < L s 
(i = 1,2,3). A gauge field U x ^ is assigned to the link pointing from point x 
to point (x + /i), where fj, — 0, 1, 2, 3 designates the four forward directions in 
space-time. 
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The order parameter (magnetization) of the system at some MC time t is the 
expectation value of the Polyakov loop 

m(<)4e^w> > (9) 

which is defined as the trace of ordered products of gauge field variables 

I Lt-l 

Lx(t) = -Tr J] U Mj0 (t) . (10) 

/ x =0 

The average is taken over independent measurements, i.e. independent random 
numbers. The deconfining phase transition of this model is of second order. 
The critical point for L t — 4 was determined with 4/g c 2 = 2.2989(1) for infinite 
large space extensions [14]. 



4 The HMC algorithm 



Let us briefly sketch the HMC algorithm. In ordinary Metropolis or heatbath 
updating algorithms the new configuration is generated by sweeping over the 
whole system and changing locally the field variables q. In case of the HMC 
algorithm one uses molecular dynamics to generate the new configurations. 
One starts by introducing an additional, fictitious so-called molecular dynam- 
ics time t' and corresponding momenta p. The initial conjugate momenta p^ 
are generated from a Gaussian distribution of unit variance and zero mean. 
The fictitious time evolution of the fields and the momenta is now given by 
the following set of coupled first-order differential equations: 

P< = — S— > Qi=Pi, (11) 



where the Hamiltonian is given by TC = Yl,iVi ! /2 + S[q\. The time derivates 
are to be understood with respect to the fictitious time if . The numerical 
integration of Eqs. (11) is performed by using a discretized version. In practice 
one uses a leap-frog integration scheme, using Nmd integration steps of size 
At' in order to integrate from fictitious time to some value t' = Nub At'. 
The endpoint of these trajectories are considered as a trial new configuration, 
which is accepted or rejected according to the general Metropolis acceptance 
probability. The HMC algorithm is exactQ, i.e. systematic errors arising from 

1 The lack of reversibility coming from round-off errors in the numerical integration 
are discussed in Ref. [16]. 
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finite time steps in the molecular dynamics are avoided by the accept /reject 
step. The algorithm is also ergodic due to the stochastic update of the initial 
momenta and fulfills the detailed balance condition, because of the reversibility 
of the leap-frog integration. 

For SU(2) lattice gauge theory the equations of motion are 

4 • 

iH XjfJi = — - (U X)tl V X)tl — h.c.) , U X>IJ , = iH x ^U x ^ , (12) 



where H X)/1 is the momentum conjugate to the field U XifM and takes the values 
in su(2), the Lie algebra of SU(2). V X)IX denotes the staples around the link 
U X fl , i.e. the incomplete plaquettes that arise in the differentiation^]. The 
classical trajectories are computed using the leap-frog scheme, which consists 
of a sequence of intermediate points (j = 0, . . . , N MD — 1) of the following form 



At' • 

H Xjlx (At'/2) = H x JO) + —H x ^(0) , 



U Xill + At') = exp (i At'H^ (j At' + ^ j j U x Jj At') 



At'\ ( Af'" 

H x>li ( j At' + — j = H,^ [j At' - — ) + At'H^U At') 



( At'\ At' ■ 
H x ,,(t') = H X4L It' -—)+ —H x At') ■ (13) 

The scheme is exact up to 0(At' 2 ). In order to generate the desired Boltz- 
mann distribution and to account for the discretization errors, the new con- 
figuration is only accepted with probability P = min{l, exp(— A7Y)}, where 
AH = H(U(t), H{t')) - H(U(0), H(0)) and the Hamiltonian is given by 

n[U,H} = W^(H x /) + S[U]. (14) 



5 Numerical results 



We perform simulations with L t = 4 and L s = 8, 16 and 24 at the critical 
point A/g c 2 = 2.2989 and in the neighbourhood^ at A/g 2 = 2.2689 and 2.3289. 

2 Details of the HMC algorithm for gauge theory can be found in Ref. [17]. 

3 For simulations in equilibrium these values would be considered to be far outside 
the critical region, especially for L s = 16 and 24. However, these values are close 
enough to the critical point to calculate the derivative of the magnetization with 
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Starting from the ordered initial state, i.e. all link variables U x ^ are set to 
unity, we measure the magnetization M and the cumulant U as a function of 
the MC time t. The system is updated with the HMC algorithm, where we fixed 
the length of the trajectory to t' = 0.32. A unit in the MC time t is defined as 
one global Metropolis step. Simulations are performed up to t = 400 global MC 
steps. The average is taken over 0(1000) samples for L s = 8, 0(100) samples 
for L s = 16 and O(10) samples for L s = 24. Statistical errors are calculated by 
dividing the data into different subsamples. Systematic errors are estimated 
by the results of different system sizes and different time intervals, i.e. we 
examined the dependency of the critical exponents from the fitted interval 
t = [tmin,£max] and the space direction L s . The quoted error is a sum of the 
statistical and systematic error. 

In Fig. 1 we plot the time evolution of the magnetization at the critical point 
for different system sizes on a double logarithmic scale. Statistical errors are 
of the order of the distance between the curves. For L s = 16 a trajectory of 
the HMC algorithm consists of 80 steps with At' = 0.004. To get comparable 
results for the other lattices (L s = 8, 24), one can not perform simulations with 
the same parameters (Nmd, A*')0- The reason is the global accept/reject step. 
The difference of the final and initial values of the Hamiltonian increases for 
larger lattices, so that using the same parameters would result in a lower 
acceptance rate. Therefore, we scale the step size with At' ~ L^ 1 to get 
comparable acceptance rates. Since the trajectory length t' is constant, we 
have to scale ^md ~ L s so that the CPU time scales also with L s . In all cases 
(L s = 8, 16, 24) we get an acceptance rate of approximately 99%. 

However, also if we scale the parameters in the way described there are large 
deviations for the different system sizes at small times. This might arise from 
statistical effects coming from the global accept/reject step. Especially for the 
largest lattice we have a small statistics and the effect of a global reject step at 
small times is high. Also, the acceptance rate (which is calculated by averaging 
over the 400 MC steps) for early times is relatively small. However, the results 
show that finite size effects coming from the finite space dimension are small 
up to t = 400, i.e. the influence of finite lattice size in Eq. (4) do not show 
up in the time interval useclQ. Therefore, L s = 16 is large enough to avoid 
systematic errors from too small lattices. 

Obviously, the magnetization which is shown in Fig. 1 can be described by 
a power law behaviour if we leave out the data for small times up to t mic . 
Thus, the time interval shown is completly within the short-time regime. The 

respect to r in the short-time region. 

4 In case of the local Metropolis algorithm simulations of different system sizes can 
be simply performed with the same parameters. 

5 A system starts to show finite size effects after a time scale proportional L z , which 
leads to deviations from the power law behaviour. 
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microscopic time scale during which the non-universal behaviour is swept away 
is about 40. The slope of the curve yields the value of the exponent C\ = /3/vz. 
The results of the different space dimensions, subsamples and time intervals 
lead to ci = 0.248(5). 

The critical exponent c\ should be independent of the step size At' and number 
of steps Nmd of the trajectory. However, a change of these parameters can 
lead to a different microscopic time scale t mic and a change of the statistical 
errors. To examine the influence of the parameters At' and N MD we study 
the short-time critical dynamics of the magnetization with At' = 0.002 and 
Nmd = 160 for L s = 16. The results are compared with the measurements 
using At' = 0.004 and Nmd = 80 and are shown in Fig. 2. The acceptance rate 
changes from 98.7% for N MD = 80 to 99.7% for N MD = 160, while the CPU 
time increases by a factor of two. Although the changes for the acceptance 
rate are small, there are large deviations between both measurements at small 
times. Figure 2 shows that the microscopic time scale t m i c decreases for smaller 
At'. The difference between both simulations at larger times are negligible and 
the values for the critical exponent c\ coincide within statistical errors. 

To extract the critical exponent Qi, we measure the magnetization as a func- 
tion of time also below and above the critical point. The simulations are per- 
formed with L s = 16, Nmd = 80 and At' = 0.004. The results are visualized 
on log-log scale in Fig. 3. These data are used to calculate the logarithmic 
derivative of the magnetization with respect to g c . This is shown in Fig. 4. 
The slope provides en = 0.83(3), where the error of c u is dominated by sys- 
tematic effects. In principle, one can also use the simulations to estimate the 
critical point. This is done by searching the best power law behaviour of M(t) 
between the two values 4/gi 2 = 2.2689 and 4/#2 2 = 2.3289 as described in Ref. 
[6]. Namely, the best straight-line fit to curves obtained by quadratic interpo- 
lation for g 1 > g > g 2 is sought. However, our statistic is not good enough so 
that we do not perform this analysis. 

The final step is to determine the critical exponent z. This is done by mea- 
suring the cumulant U(t) and extracting the exponent d/z. Our simulations 
yield z = 2.0(1). Thus we get v = 0.60(5) and (3 = 0.30(2). These results are 
in agreement (within statistical errors) with those obtained from simulations 
in equilibrium [14], which are v = 0.630(11) and (3 = 0.328(6). Also, the data 
coincide with those of the three-dimensional Ising model [7] . Performing sim- 
ilar measurements (but using the Metropolis algorithm) one gets in this case 
v = 0.6327(20), = 0.3273(17), and z = 2.042(6). 

The advantage of the short-time dynamic approach compared to simulations 
in equilibrium is that is free of critical slowing down since the spatial correla- 
tion length is small within the time regime, even at or near the critical point. 
Thus in case of local algorithms (Metropolis, heatbath) the CPU time to get 
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comparable results for different systems is independent of the system size, 
while simulations in the long-time regime suffer from critical slowing down. 
However, if we use the HMC algorithm increasing the system size leads also 
for simulations in the short-time regime to an increase of the CPU time since 
we have to scale the step size At'. Therefore, determing the static critical 
exponents can be done easier if we use local algorithms (and short-time dy- 
namics). However, if one wants to measure the dynamic critical exponent z or 
include fermions in the simulation one has to use the HMC algorithm. In this 
case using the dynamic relaxation is advantageous compared to conventional 
simulations in equilibrium. 



6 Summary and conclusions 

We presented comprehensive Monte Carlo simulations of the short-time critical 
dynamics for SU(2) lattice gauge theory in (3 + 1) dimensions. The dynamics 
of the system was given by the HMC algorithm. Starting from the ordered 
state, the magnetization (Polyakov loop), its derivative with respect to the 
coupling constant and the cumulant were measured at the critical point. The 
observables obey a power law behaviour after some microscopic time scale t mic 
as expected. The critical exponents fl/vz, 1/vz and d/z are determined from 
these time dependencies. The results support a universal short-time scaling 
behaviour for SU(2) gauge theory in (3 + 1) dimensions. The values for the 
static exponents f3 and v agree within statistical errors with those measured 
in equilibrium and with those of the three-dimensional Ising model. Thus the 
(3 + l)-dimensional SU(2) lattice gauge theory and the Ising model in three 
dimensions are in the same universality class. The dynamic critical exponent 
for the HMC algorithm was determined with z = 2.0(1). The work could 
extended to SU(2) lattice gauge theory with dynamical fermions. 
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Fig. 1. Time evolution of the magnetization at the critical point starting from the 
ordered state for L s = 8, 16 and 24. The length of a trajectory of a HMC step was 
t' = 0.32 with Nmd = 80 intermediate points. 
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Fig. 2. Magnetization as a function of time for L s = 16 and t' = 0.32 at the critical 
point. A trajectory consists of Nmd = 80 steps in the first case (full line) and 
Nmd = 160 steps in the second case (dashed line). 
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Fig. 3. Time evolution of the magnetization for three values of the coupling constant 
with L s = 16 and JVmd = 80. 
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Fig. 4. Logarithmic derivative of the magnetization with respect to r taken at g c , 
obtained from the curves shown in Fig 3. 
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